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We study a simple exactly solvable 2D model describing the interaction of a localized 
particle with an impurity. The localization potential V(x) = —aS(x) causes the particle 
to be trapped in the y-axis, and the 'impurity' is modeled by a straight impenetrable edge 
extending along the positive x-axis. We show that the problem described can be treated as 
the Sommerfeld diffraction from an infinite edge. We use the model to present qualitative 
arguments on the nature of interaction of polarized light travelling along an optical fiber 
with external impurities. 



I. INTRODUCTION 



There is a class of problems characterized by an example of polarized light propagating along 
an optical fiber (or a wave guide). It is interesting to see how the interaction with depolarizing 
impurities arises in such a problem and what effect the impurities have on the propagation of 
polarized light. The polarized light is described by the scalar Helmholtz equation. The optical 
fiber is a waveguide with appropriate boundary conditions imposed. The depolarizing interaction 
with impurities on the walls of the waveguide can be modeled by a photon spin-flip interaction. 
Since this problem is too complex to obtain an exact analytic solution, we will consider a much 
simplified model which, we believe, captures the essense of the problem. The Schrodinger equation 
for a particle of unit mass m = 1/2 is the same as the Helmholtz equation. We will model the 
waveguide by an attractive, localizing, potential V(x, y) and the impurity by an impenetrable 



barrier which we will take to be a semi-infinite edge located along the positive x-axis 



5(]. A simple 



argument based on the first order perturbation theory (the Born approximation) applied to a 
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bound state V'o i n a potential V(x, y) suggests that an "impurity" located on the x-axis at x = a, 
and modeled by a potential Vi mp (x,0) = X6(x — a)5(y), modifies the bound state wave function 
by a correction which is quadratic in the value of this wave function at the point where the 
"impurity" is located. What this means is that in the case when the bound state wave function 
decays exponentially with x as \ipo(x,y)\ ~ e~ a \ x \ the correction to the wave function caused by 
an "impurity" is of the subleading order, {x,y)\ ~ e -2aa e -a\x\ _ j n sense the bound state 
'exponential tails' are rather rigid because the "impurity" changes the asymptotic behavior of the 
bound state by an exponentially negligible correction to an overall normalization factor. 

II. THE MODEL 

Consider a particle moving in the Euclidean plane. Along the positive x-axis we place an 



^article will be 
4J. We are 



impenetrable barrier, a straight semi-infinite edge. A free quantum mechanical 
scattered by a straight edge in a manner first described by Sommerfeld [l|, [2, y, 4\. 
analyzing a situation where the particle is trapped in an attractive potential V(x,y) = —a5(x) 
that allows it to move freely only in the y direction. The Hamiltonian for this particle is the sum 
of two one-dimensional Hamiltonians H = H\ + H2 , H\ = —d% — a8(x), H2 = —dy , where a > 0. 
Therefore, the total Hamiltonian is 



H = -(d 2 x + dl) -a8(x). (1) 
In the simplest possible case that we deal with here, H\ has only one bound state, 

Mx) = e" aN • (2) 

In addition, there are scattering states, 

^( x ) = d vx + Ae~ ipx , x>0, tl>(x) = Be ipx ,x < 0, (3) 
where the transmission and reflection coefficients are 



A = ^-,B = ^. (4) 

p — ia p — ia 
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A simple pole at p = ia in the S-matrix corresponds to a bound state ([2]). 

The stationary solution of the Schrodinger equation describing a bound state of H and corre- 
sponding to the energy 



E 



a 



(5) 



is 



Mx,y) = e lky ~ alxl . (6) 

The problem we want to solve is the problem of scattering on the straight edge x > a, y = 0. We 
assume that for x > a, y = the wave function or its normal derivative must vanish, thus either 
i^(x, y) = or d y ip(x, y) = 0. These conditions are equivalent to the Dirichlet or Neumann boundary 
conditions, correspondingly. The idea that underlies the solution is to explore the Helmholtz 
equation 

(A + fc 2 )^ = 0, (7) 
in the complex x-y space (complex x and y planes), where 

A = d 2 x + d 2 y , (8) 

and x and y are independent complex variables. Elementary solutions that represent plane waves 
exp(ikix + ik2y) are holomorphic everywhere in the x and y complex planes. Since the Helmholtz 
equation is linear, the superposition principle can be used to obtain more general solutions from 
the basic ones. In particular, the superposition principle can be applied to obtain bound state 
solutions. 
Now, 

ipo = A exp(—ikr cos(ip — f3o)), (9) 

[ A(f3)exp(-ikrcos{tp- f3))d[3. (10) 
Jc 

To get ipo from the last formula we need to choose: 
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A(f3) = exp(i/3)/27r(exp(i/3) - exp(i/3 )) (11) 

and so 

, = J_ f exp(i7/2)exp(ifcrcos7) 

V 4vr 7c (exp(i 7 /2) - exp(-i X /2)) 7 ' 1 j 

where X = (^ — A))/2- For the scattering states arising in the scattering off the straight edge, 
A) = 7r/2. For the bound state propagating along the y-axis only, we have /?o = vr/2 ± zA, where 
A = X(k, a). 

The cut starting at x = a and running all the way in the positive direction of x-axis suggests 
introducing a double covering of a complex z-plane, where z = y + i(x — a). To this end, we 
introduce a new complex variable w = £ + irj defined by 

z = w 2 . (13) 

We can now write 

x — a = r cos (p = —r sin(np — vr/2), (14) 

y = r sin (p = r cos(<£> — tt/2), (15) 



and identify £ and 77 as 



£ = v 9 C os-,?7 = -v / rsm-, (16) 



where \ = <p — tt/2. We observe that both ip = and 2tt correspond to y = and x > a 
(as x = r + a), but the first case translates into £ = 77 = \/r/2, whereas the second one into 
£ = 77 = —yJr/2. For y> = 7r we obtain that £ = —r\ = y/r/2. This corresponds to y = but with 
a; < a. In these new variables the Laplacian reads 

A = 4(g2 1 +T?2) ( g l + ^ 2 )- ( 17 ) 

We will seek the solution to the scattering problem described by the Helmholtz equation ([7]) in the 
form of -0 = ip\ + ip2 1 where 

01 = exp(-iky)V(£,ri), (18) 



02 = exp(iky)U(£,ri). 



(19) 
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Plugging tpi into ([7]) yields 

"(3f + d 2 ) - Aik^ - rjd v )] V{i, rj) = 0, (20) 
and suggests that one should take V(£, rf) = V(£) and f/(£, 77) = U{rj). As a result, one obtains 

V(0 = [* dreM^kr 2 ) = F(£) = C(£) + iS(0, (21) 

J — oo 

where C(£) and S'(^) are Fresnel integrals. The expression for U(jj) is identical to that of V(£) 
except that it is in terms of rj. Consequently, 

VK£, rj) = exp(-i%)F(0 - exp{iky)F( V ), (22) 

where, as seen from (|13[) . 

y = - rj 2 , x = 2£ri + a. (23) 



For y = 0, we have 



i>(x,0) = ~ F{V) (24) 



which equals for x > a and (F(£) — F(—^)) otherwise, as should be apparent from our earlier 
discussion of the relationship between £ and rj in the case y = 0. 

To conclude this part, we have demonstrated that the solution to the pure scattering problem 
under study in the entire x-y plane is given by 



ip(x,y) = C 



r 2 P i 

exp(-iky) / drexp(2z/cr ) — exp(ifcy) / drexp(2ikr 

J — oo J —oo 



(25) 



where £ and r\ are related to x and y via ()23[) and Co is a constant. 

Let us now consider the case of the bound state and assume for simplicity that a = 0. Writing 
ipo(x,y) given by Q as ipo(x,y) = exp(iS), we obtain 

iS = —iky — a\x\ = ikr sin ip — aer cos ip 

= r ({k + as) exp(iip) — (k — ate) exp(—i(p)) /2 

= inr cos((^ — tt/2 — iX) = inr cos(c^ — xo) = lKr cos X> (26) 

where we used the substitutions A; + ae = Kexp A and k — ae = Kexp(-X). Moreover, e = e{x) = 
x/\x\ is the signum function and xo = 7r /2 + £A. The last sequence of formulas is valid for the 
positive values of energy. The expression for iS that corresponds to negative energies can be easily 



6 



obtained by analytical continuation. By analogy to the free scattering case, one can introduce 
new complex variables £ and ry related to the radial and angular coordinates r and x i n the same 
manner as in (|16h except that now 

^ = V2( C ° S ^ + Sm ^)' (27) 

ff ( (p — iX <p — iX\ 

77 = V 2 (v C ° S — ^ — ^ — J • (28) 

One can show that for <p equal and 2tt, £(93) = r/*(ip), where the asterisk denotes the complex 
conjugation. Pursuing further the analogy to the free scattering, it is straightforward to write the 
solution for the case under consideration, 

tp = exp(— a\x\) (exp(— iky)F(^) — ex.p(iky)F(r]*)) . (29) 

This is an exact solution to the scattering problem on the half-infinite edge impurity which is valid 
both for real and imaginary k. The case of imaginary k, or negative energy, is of interest because 
it corresponds to the case of a single bound state of the attrractive Dirac-delta potential in the 
transverse direction x. In this case we find an agreement with the simple perturbation theory 
argument presented in the introduction. Our investigation can be extended in many directions. In 
particular, it seems possible to extend the method presented here to localized wave packets. 
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